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We study the evolution of massless scalar waves prop- 
agating on spherically symmetric spacetimes with a non- 
zero cosmological constant. Considering test fields on both 
Schwarzschild-de Sitter and Reissner-Nordstrom-de Sitter 
backgrounds, we demonstrate the existence of exponentially 
decaying tails at late times. Interestingly the I = mode 
asymptotes to a non-zero value, contrasting the asymptoti- 
cally flat situation. We also compare these results, for £ — 0, 
with a numerical integration of the Einstein-scalar field equa- 
tions, finding good agreement between the two. Finally, the 
significance of these results to the study of the Cauchy horizon 
stability in black hole-de Sitter spacetimes is discussed. 



I. INTRODUCTION 

Perturbative studies of relativistic, spherical collapse 
have elucidated dynamical features of gravitational col- 
lapse important to an understanding of black hole forma- 
tion and the subsequent relaxation to a stationary state 
(see for example jlj). Indeed, quasi- normal ringing could 
provide direct evidence of the existence of black holes 
if observed §]. At late times, quasi-normal oscillations 
are swamped by the radiative tail of the gravitational 
collapse |2|] . This tail radiation is of interest in its own 
right since it originates from the non-trivial propagation 
of the field perturbations on the background spacetime — 
zero rest-mass fields do not necessarily propagate along 
the light-cone in a curved spacetime. 

The first authoritative study of nearly spherical col- 
lapse, exhibiting radiative tails, was performed by Price 
j|. Studying the behaviour of a massless scalar field 
propagating on a fixed Schwarzschild background, he 
showed that the field dies off with the now familiar power- 
law tail t-( 2e + p + 1 ) a t late times, where I is the multipole 
order of the field, and P = 1 if the field is initially static 
and P = 2 otherwise. Furthermore, Price showed that 
the perturbations of any zero rest-mass, integer-spin field 
are governed by a wave equation with the same qualita- 
tive form as that governing the scalar field. This suggests 
that the results for the scalar field should apply equally 
well to the radiatable multipoles of both the electromag- 
netic and gravitational fields. Similar results for a mass- 



less scalar field propagating on a Reissner-Nordstrom 
background have been obtained by Bicak || . No such an- 
alytic result has yet been obtained for the case of a black 
hole with angular momentum, though Krivan et al. || 
have recently performed numerical work which suggests 
the power law tail holds independently of the angular 
momentum of the black hole. 

While test-field calculations arc extremely compelling, 
it is natural to ask to what extent linear analyses are 
representative of dynamical gravitational collapse. If ci- 
ther quasi-normal ringing, or radiative tails, should be 
absent in non-linear collapse one might view results of 
linear analyses with scepticism. Advances in numeri- 
cal relativity make it possible to address this issue in 
the spherically symmetric context. Gomez and Winicour 
H studied the non-linear evolution of a self-gravitating, 
spherically symmetric, massless scalar field concluding 
that the scalar monopole moment decayed exponentially 
rather than with the power-law predicted by the linear 
analyses. More recently, Gundlach, Price and Pullin [QJ|] 
re-examined this problem. They were able to show that 
the frequencies of quasi-normal oscillations, and decay 
rates of the power-law tails, found in the numerical so- 
lutions, arc in good agreement with the predictions of 
perturbation theory (though one must go to sufficiently 
late times in order to see tail effects - this would, in part, 
explain the null result of Gomez and Winicour) . 

The presence, and slow decay, of wave tails at late 
times is a key ingredient leading to the instability of 
Cauchy horizons inside charged and rotating black holes. 
For black holes in asymptotically flat spacetime, the in- 
verse power-law decay of perturbing fields at the event 
horizon has been used to provide initial data in linear |{J 
and non-linear [p~0|— p~3[ studies of the black hole interior. 
In particular, the form of the wave tail is largely respon- 
sible for the weakness of the mass-inflation singularity 
inside charged black holes and is believed to have 
similar consequences inside rotating holes ]l2| ] . 

In contrast to the asymptotically flat case, perturba- 
tive |l4|-|l6|| and non-linear [|l7| studies indicate that the 
inner Cauchy horizon of charged and rotating black holes 
immersed in asymptotically de Sitter space can be stable. 
That stability persists for a finite volume of the parame- 
ter space suggests that these spacetimes violate the letter 
(if not the spirit) of the cosmic censorship hypothesis Jl8| . 
The nature of the radiative tail of perturbations at late 
times plays a major role in these analyses. Plausible ar- 
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guments suggest that an exponential decay of the tails 
replaces the power law behaviour observed in asymptot- 
ically flat space |l5||l6|], however, no detailed analysis of 
the evolution of wave tails in asymptotically de Sitter- 
black hole spacetimes exists. The present work reports 
on such a study for non-rotating black holes. Our pri- 
mary motivation has been to obtain correct boundary 
conditions on the radiation at the event horizon for use 
in a numerical study of the internal structure of charged 
black holes in de Sitter space ||l6|| . 

The paper is organized as follows: In section II we 
consider the propagation of massless, minimally coupled 
scalar fields on the Schwarzschild-de Sitter and Reissner- 
Nordstrom-de Sitter black-hole backgrounds. We derive 
the equation governing the scalar test-field, and numer- 
ically integrate it. Two independent numerical codes 
were used throughout the linear analysis; a null evolu- 
tion scheme following that of Gundlach et al. (7), and 
a Cauchy evolution scheme similar to that described by 
Krivan et al. We found complete agreement between 
them. Our results show that, except for the I = mode, 
the field falls off exponentially at the cosmological and 
black-hole event horizons, and at future timelike infinity. 
The rate of decay depends upon the surface gravity, K\, of 
the cosmological horizon [see Eq. (Q)], and the multipole 
order of the field. In particular, as t — > 00 

4>~e- tKlt , £>0, (1) 

where t is defined by Eq. (^). For I — 0, the field ap- 
proaches a constant value at late times rather than de- 
caying. A suggestion of this behaviour can be found in 
the analysis of Chambers and Moss j[6| and, as argued 
there, is similar to the situation within the black hole 
interior. In section III, we study the non-linear evolution 
of a spherically symmetric, self-gravitating scalar field by 
numerically integrating the coupled Einstein-scalar field 
equations. Confining attention to spherical symmetry 
implies we gain information solely about the I = mode 
of the field. We find, in accord with our linear analysis, 
that the field approaches a constant value at the cosmo- 
logical event horizon, the black hole event horizon and 
future timelike infinity, demonstrating that the results 
of our linear analysis are indicative of the full theory. 
Further to this, we inspect the behaviour of the field's 
stress-energy, showing that 

<f>,r ~ e~ 2ku , (2) 

with k — K± to within about 12%. In section IV we make 
some final comments about the implications of our re- 
sults for Cauchy horizon stability in black-hole-de Sitter 
spacetimes and the related issue of cosmic censorship. 

II. A LINEAR ANALYSIS 

In this section, we consider a massless, minimally 
coupled, scalar field propagating on a fixed Reissner- 



Nordstrom-de Sitter background. Since our considera- 
tions are limited to the black hole exterior, it is clear 
that the Schwarzschild-de Sitter case can always be ob- 
tained by setting the charge q = (when this fails, we 
explicitly include the corresponding formulae). 

A. The equations 

The generalization of the Reissner-Nordstrom metric 
to include a cosmological constant was first given by 
Carter |J] as 

ds 2 = -f{r)dt 2 + fir^dr 2 + r 2 (d9 2 + sin 2 6d<p 2 ) , (3) 
where 

/(r) _ 1 _^ + 4_4 . * » >0 . (4) 
r r z or A 

In Eq.(Q), M denotes the mass of the black hole, q its 
charge and A is the cosmological constant. If q ^ there 
are three horizons located at the roots of f(r) = 0; an 
inner (Cauchy) horizon at r = r%, a black hole event 
horizon at r = r-i and a cosmological horizon located at 
r = 7"i such that r\ > T2 > r% (see Fig. |J). The fourth 
root, r4, is negative and thus non-physical. (When q = 
there are only two horizons (see Fig. |l|) the black hole 
event horizon, T2, and the cosmological event horizon, 
n, while the third root, r$ is negative.) 



r=0 j+ 
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FIG. 1. A conformal diagram representing the 
Schwarzschild-de Sitter black hole spacetime. Shown are the 
cosmological horizon at r = r%, the black hole event horizon 
at r = T2 and the singularity (wavy line) located at r = 0. 
Also shown are the locations of past and future null infinity 
J~ and J + respectively. 
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It is convenient to introduce a "tortoise" radial coordinate^here^^^/^Tdjecwmpbstato 

multipole pieces, i.e. <j> — J2 ^( u i v )Ytm{0, 4>)r~ l . The 
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where the arbitrary constant of integration has been 
set to zero. 




FIG. 2. A conformal diagram showing the Reiss- 
ner-Nordstrom-de Sitter black hole spacetime. Shown are the 
cosmological horizon at r = ri, the black hole event horizon 
at r — T2, the inner (Cauchy) horizon at r = r$ and the sin- 
gularity (wavy line) located at r — 0. Past and future null 
infinity are indicated by J~ and J + respectively. 



We define 



df(r) 



dr 



(6) 



where are the roots of f(r) = 0. When the root corre- 
sponds to a physical horizon in the spacetime, Ki is the 
surface gravity of that horizon p(J . Finally, we introduce 
a pair of null coordinates on the spacetimes, the advanced 
time v = t + r„ and the retarded time u — t — r*, in terms 
of which the interval (^) reduces to, 



ds 1 



-f(r)dudv + r 2 (d9 2 



(7) 



The definition we have adopted means that the future 
cosmological horizon r = r\ is located at v = oo, and 
the future black- hole event horizon r = r 2 is at u = oo. 
In terms of the null coordinates (u,v), the scalar wave 
equation, □</> = 0, becomes 
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(8) 



is highly localized near to r* =0, falling off exponentially 
in at both r = t\ and r = r 2 . The form of the potential 
for ^ = 0,1 and I — 2 is shown in Fig. |[ 




FIG. 3. The 'effective' potential when q = 0.5, M = 1, 
A = 1(T 4 and £ = (solid), i= \ (dashed) and i = 2 (dot- 
ted). We have multiplied by |r — r2\/£\ in order to accentuate 
the nature of the potential when £ = — the usual potential 
barrier is followed by a potential well, a feature not evident 
when £ > 0. The potential tends to zero exponentially quickly 
as r, — » ±oo. 



B. The results 

It is straightforward to integrate Eq. (||) on a null grid 
using the methods described by Gundlach et al. jjj . Fur- 
ther details can be found in their article. 

Since tail effects are independent of the initial data 
used, we chose to represent a generic initial perturbation 
by a Gaussian pulse onu = 



^{u = 0, v) = Aexp{-(w - vi) 2 /a 2 } 



(10) 



(The amplitude A is irrelevant since Eq. (g) is linear. The 
data used to produce the figures had center v\ — 50.0, 
and width a = 3.0.) The field is constant on v = 0, 
^(v = 0) = ^(u — Q,v = 0). We have set the mass of 
the hole, M, equal to unity throughout; this corresponds 
to the freedom to rescale the coordinates by an overall 
length scale. Investigations have shown that the results 
are qualitatively similar for all A > 0, so we fix A = 10~ 4 
from here on. We discuss the behaviour of the field, <fi — 
'F/r, in three regions: (a) timelike infinity — approached 
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on surfaces of constant r. (b) The cosmological horizon — 
in practice, approximated by the null surface v = i> ma x, 
the largest value of v on our grid, (c) The black-hole 
event horizon — again, approximated by the null surface 
u = Umax; the largest value of u on our grid. 

Gundlach et al. argued that the nature of the tails 
in asymptotically flat spacetimes is primarily due to the 
power-law form of the effective potential as r, — > oo. 
Moreover, Ching et al. |2lJ] have demonstrated that in- 
verse power-law tails, as seen in spherically symmetric, 
asymptotically flat spacetimes, are not generic. There- 
fore, it should not be surprising that tails, in asymptoti- 
cally de Sitter spacetimes, fall off exponentially with time 
since the effective potential is exponentially suppressed 
as r* — > ±00. What is unexpected is that the the 1 = 
modes do not decay to zero, rather a generic perturba- 
tion leads to a residual constant field at late times. (Of 
course, there is no stress-energy associated with a con- 
stant field.) This was suggested in a paper by Chambers 
and Moss |f6| where it was argued that the situation in 
Schwarzschild-de Sitter spacetime is somewhat analogous 
to scattering inside a black hole |^| , and hence a constant 
mode can be transmitted to both the black-hole and cos- 
mological event horizons. Indeed, the field at late times 
behaves like 

<f>\ l=0 ~ 0o + Mr)e~ 2Klt ■ (11) 




FIG. 4. A plot of \(f>\ versus time for q = 0.5, i = 0, 
M = 1 and A = 10 -4 . The field is shown on the cosmo- 
logical event horizon ri (short-dashed), the black hole event 
horizon r-2 (dotted) and two surfaces of constant radius (solid 
and long-dashed). Initially the quasi-normal ringing domi- 
nates, while at late times the field settles down to the same 
constant value on all four surfaces. 

Figure [| demonstrates this effect. The early time be- 
haviour of the field, </>, is dominated by the quasi-normal 
ringing, whose explicit features we are not currently in- 
terested in. At later times the field quickly approaches 
the same constant value at the cosmological event hori- 
zon, the black-hole event horizon and future timelike in- 
finity. 

One might suspect that the existence of a non-zero field 
at late times is an artifact of the initial data, or that 
setting q = might produce different results. There- 
fore Fig. H shows results for a Schwarzschild-de Sitter 
black hole with i = 0. Once again, the field settles 
down to a constant value which is independent of the 
radial position. One must contrast this to previous re- 



sults for asymptotically flat spacetimes where the field 
always approaches zero. Monitoring the behaviour of 
(j) t t during the evolution (this is a matter of formality 
in the Cauchy code, as it uses momenta pt = 4>,t for 
the evolution) has allowed us to conclude that the nu- 
merical value of the exponent in Eq. (|ll]) is indeed K\ to 
within 10%: Taking the time derivative of Eq. (ftl|) yields 
In 1 >t | = ln(2Ki(/)i(7-)) — 2ki£, thus the fitted values for 
Ki can be compared with the ones obtained from Eq. (||) . 

In general, the final field value, <fio, is a function of 
the black hole parameters (M, q, A). We investigated the 
dependance of </>o on the cosmological constant A. As 
shown in Fig. |6] for the case q = 0,M = 1, our results 
indicate that <fio ~ A 6 . The linear correlation coefficient 
Ciin of In |0o| versus In A is given by ci; n = 0.98. 
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FIG. 5. A plot of \ij>\ versus time when q = 0.0, £ = 0, 
M = 1 and A = I0" 4 . The field is exhibited on the same 
surfaces as the q = 0.5 cases. Once again the asymptotically 
constant field is evident. 
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FIG. 6. A plot of 1 0o | versus A for q = 0.0, I = 0, M = 1. 
The linear correlation coefficient cvj n of In |0o| versus In A is 
given by crj n = 0.98. The dashed line represents the least 
square fit to |0 O | = b A a with a = 0.59, b = 0.09. 

For £ > the picture is different, as seen from Fig. |?] 
and Fig. ||. The early time behaviour of the field is still 
dominated by complicated quasi-normal ringing, but at 
late times a definite exponential fall-off is manifest. In 
particular, the late time wave tails are well approximated 
by ~ exp(— init). In a series of separate evolutions 
we found that for sufficiently small values of A, a regime 
of power law decay followed the quasi-normal ringing, 
though at late times the exponential decay described 
above always dominates. 



FIG. 7. A plot of |0| versus time for q = 0.5, 1 = 1, M = 1 
and A = 10 -4 . The field along is shown on the cosmological 
event horizon n (short-dashed), the black hole event horizon 
T2 (dotted) and two surfaces of constant radius (solid and 
long-dashed). The field falls off as exp( — kt) at late times, 
with k ~ Ki to within about 2%. Note that the ordinate scale 
is logarithmic. 
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FIG. 8. A plot of |0| versus time for q = 0.5, 1 = 2, 
M — 1 and A = 10~ 4 . The held is shown on the cosmological 
event horizon n (short-dashed), the black hole event hori- 
zon T2 (dotted) and two surfaces of constant radius (solid and 
long-dashed) . At early times quasi-normal ringing completely 
dominates, but eventually the field falls off as exp (— 2kt). The 
value of k ~ ki is again accurate to approximately 2%. Note 
that the ordinate scale is logarithmic. 
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III. NON-LINEAR ANALYSIS 

Given the somewhat unusual behaviour of the £ = 
modes elucidated by the test field analysis, it seems 
necessary to examine the non-linear evolution of a self- 
gravitating massless, scalar field in the presence of a cos- 
mological constant. This situation is described by the 
coupled Einstein-scalar field equations 

G a p — BnT a p — 3 Q /3A , (12) 

where 

T a/ 3 = 4>,<x4>,p - (l/2)g Q/3 (0 i7 0' 7 ) , (13) 

is the stress-energy of the scalar field which satisfies D<j> = 
0. Restricting attention to spherical symmetry (thus we 
only gain information about the £ = mode) the line 
element can be written as 

ds 2 = -ggdu 2 - 2gdudr + r 2 (d8 2 + sin 2 dd<\> 2 ) , (14) 

where g = g(u, r), g — g(u, r) and u is the retarded time. 
The coordinates have been normalized so that u is the 



5 



proper time at the origin, thus <?(0, r) = g{0,r) = 1. It 
is customary p2| to introduce two new fields h(u, r) and 
h(u, r) defined by 



= h = 



h dr . 



(15) 



In terms of these variables the field equations (|12|) be- 
come 



(lng), r = Anr^ih - h) 2 
NXr = - Ar 2 ) , 
(rh)^ r = h , 



and the wave equation is 
,« 2 ,r 2r 



".9(1 - Ar 2 ) 



(16) 
(17) 
(18) 



(19) 



A well established numerical algorithm exists to integrate 
these equations H p3[]24|l . We refer the reader to these 
articles for details. 



A. The results 

Initial data for the non-linear equations is supplied on 
an initial null cone centered on the origin. We considered 
a Gaussian pulse on u = 0, with an amplitude </>o( r / fo) 2 > 
width er, and centered on r = fo. The code was tested 
against the exact solution in p5| . We also reproduced 
power-law tails in complete agreement with [|| , when A = 
0. 

Figure ^ shows h (or equivalently (j>) at the cosmologi- 
cal event horizon, the black hole event horizon and along 
an r = constant surface. The agreement with the linear 
analysis is remarkable. We again see an initial period 
of quasi-normal ringing which decays, leaving behind a 
constant field at late times. 
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FIG. 9. A plot of |0| versus time for A = 1CT 4 . The field is 
shown along three surfaces: the cosmological event horizon ri 
(short-dashed), the black hole event horizon T2 (long-dashed) 
and a surface of constant radius (solid). After a period of 
quasi-normal ringing the field settles down to the same con- 
stant value at all three surfaces. 



Furthermore, using ([L5|) and ( |18| ) we can write r0, r = 
(h — h) and examine the stress-energy of the field at late 
times. Figure |l^ suggests that (h — h) ~ exp[— 2k\u\ at 
late times, as suggested by the the test-field analysis. 
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FIG. 10. A plot of \<f>, r \ versus time for A = 10 -4 along the 
same three surfaces. At late times exponential decay with the 
approximate form exp (— 2ku) is evident — we have determined 
the value of k ~ K\ within about 12%. 



IV. CONCLUDING REMARKS 

This study of wave-tail evolution in black hole-de Sit- 
ter spacetimes has revealed consistent but interesting re- 
sults. The radiative tails associated with a massless, 
minimally coupled scalar field propagating on the fixed 
backgrounds of a Schwarzschild-de Sitter and Reissner- 
Nordstrom-de Sitter black hole decay exponentially at 
the cosmological horizon, the black hole event horizon 
and at future timelike infinity. That is, for all modes 
except the I = mode, it approaches a constant at late 
times. We have further explored this result by consid- 
ering the non-linear problem of a scalar field coupled to 
gravity and numerically integrating the Einstein-scalar 
field equations. These results show excellent agreement 
with the linear analysis. 

Though the introduction of a non-zero cosmological 
constant into the Einstein equations may be argued to be 
somewhat unrealistic, it is the issue of Cauchy horizon 
stability in black hole-de Sitter spacetimes ]l4|-|l6|] and 
the related issue of cosmic censorship |ll| , that motivates 
this investigation. The power-law tails found by Price |!|| 
have been used as initial data in arguments pertaining 
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to the instability of the inner (Cauchy) horizon of both 
the Reissner-Nordstrom and Kerr black hole spacetimes 
and the associated phenomenon of mass inflation. Lin- 
ear perturbation studies suggest that the Cauchy horizon 
inside black holes embedded in de Sitter space may be 
stable — fluxes of linear perturbations remain bounded at 
the Cauchy horizon. This in itself does not guarantee 
stability though. A fully non-linear analysis of the black 
hole interior is needed, for which the late-time wave tails 
at the event horizon serve as initial data. With the results 
of this paper in hand, we can now embark on a numeri- 
cal study of the Cauchy horizon in Reissner-Nordstrdm- 
de Sitter, along the lines of Brady and Smith |D|. 
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